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ABSTRACT 

Results are reported from a search for low-multiplicity neutrino bursts in the Sudbury Neu- 
trino Observatory (SNO). Such bursts could indicate detection of a nearby core-collapse supernova 
explosion. The data were taken from Phase I (November 1999 - May 2001), when the detector 
was filled with heavy water, and Phase II (July 2001 - August 2003), when NaCl was added to the 
target. The search was a blind analysis in which the potential backgrounds were estimated and 
analysis cuts were developed to eliminate such backgrounds with 90% confidence before the data 
were examined. The search maintained a greater than 50% detection probability for standard 
supernovae occuring at a distance of up to 60 kpc for Phase I and up to 70 kpc for Phase II. No 
low-multiplicity bursts were observed during the data-taking period. 

Subject headings: supernovae: general 

1 



1 Department of Physics, University of Alberta, Edmon- 
ton, Alberta, T6G 2R3, Canada 

2 Department of Physics and Astronomy, University of 
British Columbia, Vancouver, BC V6T 1Z1, Canada 

3 Chemistry Department, Brookhaven National Labo- 
ratory, Upton, NY 11973-5000 

4 Ottawa-Carleton Institute for Physics, Department of 
Physics, Carleton University, Ottawa, Ontario K1S 5B6, 
Canada 

5 Physics Department, University of Guelph, Guelph, 
Ontario NIG 2W1, Canada 

6 Department of Physics and Astronomy, Laurentian 
University, Sudbury, Ontario P3E 2C6, Canada 

7 Institute for Nuclear and Particle Astrophysics and 
Nuclear Science Division, Lawrence Berkeley National Lab- 
oratory, Berkeley, CA 94720 

8 Los Alamos National Laboratory, Los Alamos, NM 
87545 

9 Laboratorio de Instrumentacao e Fi'sica Experimental 
de Parti'culas, Av. Elias Garcia 14, 1°, Lisboa, Portugal 

10 Department of Physics and Astronomy, Louisiana 
State University, Baton Rouge, LA 70803 

11 Laboratory for Nuclear Science, Massachusetts Insti- 
tute of Technology, Cambridge, MA 02139 

12 Department of Physics, University of Oxford, Denys 
Wilkinson Building, Keble Road, Oxford 0X1 3RH, UK 

13 Department of Physics and Astronomy, University of 
Pennsylvania, Philadelphia, PA 19104-6396 

14 Department of Physics, Queen's University, Kingston, 
Ontario K7L 3N6, Canada 

15 Rutherford Appleton Laboratory, Chilton, Didcot 
OX11 0QX, UK 

16 Center for Experimental Nuclear Physics and Astro- 
physics, and Department of Physics, University of Wash- 
ington, Seattle, WA 98195 

17 SNOLAB, Lively, ON P3Y 1N2, Canada 

18 Department of Physics, University of Texas at Austin, 
Austin, TX 78712-0264 

19 TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T 
2A3, Canada 

20 Present address: Center of Cosmology and Particle 
Astrophysics, National Taiwan University, Taiwan 

21 Present address: Center for Astrophysics and Space 
Astronomy, University of Colorado, Boulder, CO 

22 Present address: Department of Physics, University of 
Chicago, Chicago, IL 

23 Present address: Department of Physics, Hiroshima 
University, Hiroshima, Japan 

24 Present address: Sanford Laboratory at Homcstakc, 
Lead, SD 

25 Present address: Department of Physics, University of 
Liverpool, Liverpool, UK 

26 Present address: Dept. of Physics, University of Cali- 
fornia, Santa Barbara, CA 

27 Present address: Department of Physics, Queen's Uni- 



1. Introduction 

Supernova neutrinos offer unique insights into 
both the fundamental nature of neutrinos and 
the complex process of core collapse. Although 
theoretical calculations predict that a typical 
supernova releases approximately 3 x 10 53 ergs 
of gravitational binding energy, 99% of which 
is carried away by neutrinos, the only super- 
nova neutrinos ever det ected came from a single 
supernova, SN 198 7 A (|K. S. Hirata et all [l987t 



R. M. Bionta et al. 



19871 



E. N. Alcksccv et al. 



19871 ). Many open questions in supernova core- 
collapse models could be resolved with additional 
supernova neutrino data, motivating large neu- 
trino detectors to search their datasets for multi- 
ple events clustered closely in time, which could be 
considered candidate supernova neutrino events. 

One such neutrino detec tor, the Sudbury Neu - 
trino Observatory (SNO) (|j. Boger et aLl boOOt) . 
was an imaging water Cherenkov detector located 
at a depth of 5890 m of water equivalent in the 
Vale Inco., Ltd. Creighton mine near Sudbury, 
Ontario, Canada. SNO detected neutrinos using 
an ultra-pure heavy water ( 2 H20, hereafter D2O) 
target contained in a transparent spherical acrylic 
vessel 12 m in diameter. Neutrino interactions in 
the vessel produced Cherenkov light that was de- 
tected by an array of 9456 photomultiplier tubes 
(PMTs) supported by a stainless steel geodesic 
structure. Outside the geodesic structure were 
5700 tonnes of light water that were monitored by 
outward-facing PMTs in order to identify cosmic- 
ray muons. The D2O target allowed SNO to detect 
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neutrino events via three different reactions: 



v x + e 
v e + d 
v x + d 



->■ v x + e 
p + p + e~ 
■p + n + u x 



(ES) 
(CC) 
(NC) 



The neutral current (NC) reaction allowed SNO 
to detect all active neutrino flavors with equal sen- 
sitivity, while the charged current (CC) reaction is 
exclusive to electron neutrinos. The elastic scat- 
tering (ES) reaction is primarily sensitive to elec- 
tron neutrinos; other flavors can undergo ES reac- 
tions but with a smaller cross section. Addition- 
ally, SNO could see electron antineutrinos through 
inverse /3 decay on deuterium and hydrogen, as 
well as through elastic scattering: 



v f + e " 



d 
P 



n + n + 
n + e + 
-> i7 P + e: 



The first reaction can provide a triple coincidence 
between the two neutrons and the positron, but it 
has a relatively small cross section. The second re- 
action has a much higher cross section but occured 
only in SNO's 1700 tonne light water shield, which 
was between the acrylic vessel and the PMT sup- 
port structure. For this analysis we have focused 
only on the D2O region. Additional reactions on 
oxygen isotopes are also possible, but the rarity 
of 17 and 18 in both the H 2 and D 2 vol- 
umes made the event rates from these processes 
very low. 

SNO's sensitivity to all neutrino flavors and 
the comparison of the rates of the different pos- 
sible reactions provide exciting opportunities to 
distinguish various supernov a models and inves- 
tigate neutrino properties (K. Scholberd 2007t 



W. C. Haxton 1120081: IR. C. Schirato and G. M. Fulled 
20021 ). During SNO's run time, it participated in 
the SuperNova Early Warning System (SNEWS), 
which was designed to notify the astronomical 
community within minutes of the dete ction of a 
large neutrino burst (|P. Antonioli 2004 ) . In order 
to avoid triggering on a false burst during this 
real time analysis, the threshold for the number of 
events qualifying as a burst was conservatively set 
to be 30 events in less than 2 s. Although SNO did 
not observe such a large burst during its run time, 
a more thorough search is required to determine 



whether or not SNO observed any evidence of a 
supernova. 

In this paper we have searched the SNO dataset 
for low-multiplicity bursts, which we define as 
bursts of two or more events. Such bursts could 
come from supernovae in satellites of the Milky 
Way that are hidden by interstellar dust, from 
non-standard supernovae in our own galaxy with 
relatively low neutrino fluxes, or from completely 
unknown and unexpected sources of neutrinos. 
The Supcr-Kamiokande coll aboration has recentl y 
published a similar search ( M. Ikeda et al. l2007t ) . 
Our search is complementary to that of Super- 
Kamiokande, in that for much of the running pe- 
riod examined here the Super-Kamiokande detec- 
tor was not operational. Additionally SNO was 
primarily sensitive to electron neutrinos, while 
Super-Kamiokande was primarily sensitive to elec- 
tron antineutrinos. 

2. Data Analysis 
2.1. Data Set 

The data analyzed here include two phases of 
SNO's operation. Phase I ran from November 
1999 to May 2001, and the sensitive volume of 
the detector was filled only with D2O. Phase II 
ran from July 2001 to August 2003, and during 
this phase NaCl was added to the detector, in- 
creasing the sensitivity to the NC reaction through 
the consequent enhancement of neutron detection 
efficiency. Phase II began running shortly af- 
ter Super-Kamiokande's first phase of data taking 
ended, meaning that the majority of SNO Phase II 
contains no overlap with the supernova search per- 
formed by Super-Kamiokande. Because of the 
enhanced NC detection efficiency, SNO Phase II 
provides a higher sensitivity to a potential super- 
nova signal. The total livetime of Phase I was 
241.4 days, while the total livetime of Phase II 
was 388.4 days. The absolute time of an event 
was measured by a GPS system whos e accuracy 
was ~ 300 ns (|B. Aharmim et al\\200l\) . 



2.2. Search Parameters 

We estimated all of the backgrounds that could 
mimic a supernova burst signal, and we designed 
our search windows and analysis cuts to ensure 
that we were 90% confident we would not see 
a false burst. We optimized our searches by 
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adjusting the minimum burst multiplicity, the 
length of the coincidence time interval, and the en- 
ergy threshold for individual events. Our energy 
threshold sets a lower limit on the reconstructed 
total energy of the detected electrons, which may 
have significantly lower energies than the neutri- 
nos themselves. Sensitivity to a supernova signal 
is increased by selecting a long coincidence time 
interval, a low multiplicity requirement, and a low 
energy threshold, but these criteria must be bal- 
anced against the increase in the probability of a 
"background" burst, especially through accidental 
coincidences. As the energy threshold is lowered, 
the event rate goes up, and hence the acciden- 
tal coincidence rate increases, particularly for low- 
multiplicity bursts and for a long coincidence time 
interval. To ensure a 90% probability of seeing no 
false bursts within a given search window, the ex- 
pected backgrounds from all sources must sum to 
no more than 0.11 events. Our optimization ulti- 
mately led us to perform three distinct searches: 
two multiplicity two (iV eV ent — 2) searches, and 
one multiplicity three (iV even t = 3) search. 

The first of our -/V cvcnt = 2 searches used 
a short window (0.05 s) to focus on detect- 
ing neutrinos from a supernova neutronization 
burst. In the case of a failed supernova, the 
neutronization burst provides the only poten- 
tial signal because shortly after the neutron- 
ization phase, the supernova collapses into a 
blac k hole, abruptly terminating the neutrino sig- 
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2006). These unusual supernovae are of special 
interest to astronomers, and their neutrino signa- 
tures could provide interesting model constraints. 
The second multiplicity two window was of moder- 
ate length (0.2 s for Phase I and 1 s for Phase II) to 
maximize our sensitivity to a standard supernova 
event. Table [T] summarizes the search windows 
and their respective energy thresholds. 

The Phase I data set entirely overlaps with 
the running of Super-Kamiokande, while much 
of Phase II does not overlap. Consequently for 
Phase I our optimization was intended to main- 
tain some neutral current sensitivity, which Super- 
Kamiokande does not have to great extent, with 
the aim of being able to detect non-standard 
burst sources. For Phase II, however, our goal 



was to maximize overall supernova sensitivity, and 
therefore we used a fairly large (Is) search win- 
dow. The 1 s window required us to raise the 
energy threshold high enough that there is very 
little remaining neutral current signal, but it in- 
creased our overall sensitivity to standard super- 
novae bursts. 

For the multiplicity three search, where acci- 
dentals are not a major background, the window 
for each phase is 10 s, and the energy threshold 
is 4.5 MeV. Primarily because of the low energy 
threshold available for this search, it provides the 
best sensitivity to standard supernovae bursts. 

In order to simulate a "standard" supernova 
burst, we utilized supernova Mon te Carlo simula- 



tions based on the Burrows model (|A. Burrows. D. Klein, and R. Ganc 

19921 ). The average neutrino energies from 

this model are (E U<1 ) ~ 13 MeV, (EvJ ~ 

15.5 MeV, and (E„ v ,,) ~ 20 MeV, where 

rep resents any of the flavors v a , F u , v T , or 

v T (T. A. Thompson. A. Burrows, and P. A. Pintol 

2003b . 



2.3. Backgrounds 

The primary difficulty in performing a trigger- 
less burst search is the elimination of the "back- 
ground" created by false bursts. In addition to the 
accidental coincidences, a large number of corre- 
lated physics backgrounds must be estimated and 
almost entirely eliminated. SNO's sensitivity to 
neutrons makes this problem particularly difficult: 
any process that produces multiple neutrons can 
lead to an apparent burst, with the average time 
between neutron captures being roughly 50 ms in 
Phase I due to capture on deuterium and 5 ms 
in Phase II due to capture on CI. Table [2] shows 
our background estimates for each of our search 
windows before we apply any special analysis cuts 
designed to remove false bursts beyond the stan- 
dard SNO analysis cuts ( B. Aharmim et all feoiO. 
20071 ) . Most of the multiplicity three backgrounds 
are conservatively assumed to have upper limits 
corresponding to the multiplicity two estimates. 

The dominant correlated backgrounds are in- 
teractions by atmospheric neutrinos, which can 
produce neutrons without any primary large en- 
ergy deposit to tag the events. For both Phase I 
and Phase II, the atmospheric background was es- 
timated with the neutrino interaction generator 
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Table 1 

Search windows and energy thresholds 







Window Length 

(s) 


Energy Threshold 
(MeV) 


Phase I, N ovcnt 


= 2 


0.05 


5.0 


Phase I, N ovcnt 


= 2 


0.2 


6.0 


Phase I, Ncvont 


= 3 


10.0 


4.5 


Phase II, N vont 


= 2 


0.05 


6.5 


Phase II, N eV ent 


= 2 


1.0 


8.5 


Phase II, Ncvont 


= 3 


10.0 


4.5 



Note. — The above values were chosen to maximize our sensi- 
tivity while limiting our probability of seeing a false burst. The 
energy threshold sets a lower limit on the reconstructed total 
energy of the detected electrons. 



NUANCE (|D. Casperl 120021) . whose output was 
then further processed by the full SNO detector 
simulation. The simulation's atmospheric neu- 
trino energies ranged from 100 MeV to 2 TeV, and 
flavor oscillation corrections were applied. The 
systematic error in the NUANCE simulation is 
conservatively estimated to be ±20%, and it is 
dominated by uncertainties in the neutrino cross 
sections. 

Most muons traveling through the SNO detec- 
tor were tagged by outward-looking PMTs, and 
neutrons following these muons are eliminated by 
a 20 s software 'muon follower' veto. Some muons, 
however, do not have enough energy to trigger the 
outward-looking PMTs and will leak into the de- 
tector. Fortunately few of these muons are likely 
to produce multiple neutrons that could mimic a 
burst. In the entire livetime of Phase II, the ef- 
fect of the remaining leaked muons is estimated to 
cause fewer than 1.35 single- neutron events, im- 
plying a conservative upper limit of fewer than 
0.5 coincidence events for Phase II. Since Phase I 
would have even fewer coincidences, the same con- 
servative upper limit of 0.5 coincidences is as- 
sumed. 

Although the SNO detector was constructed 
from materials with low radioactivity, some ra- 
dioactive backgrounds still existed, such as 238 U, 
and spontaneous fission from residual radioactiv- 



ity can lead to false bursts due to multiple neutron 
capture. Many of the radioactive backgrounds dis- 
cussed in previous SNO analyses are not signifi- 
cant in this search because they will not produce 
bursts of events, but fission neutrons from 238 U 
can create a background burst. The amount of 
238 U in the detector was measured using an in- 
ductively coupled plasma mass spectrometer in 
September 2 003 after the addition of NaCl to the 
heavy water (|B. Aharmim et aZ.ll2004h . More 238 U 
could have been present in Phase II than in Phase I 
because the addition of NaCl prevented the use of 
SNO's reverse osmosis purification system; there- 
fore, this 238 U measurement was used as an upper 
limit for Phase I. Using the standard probability 
of a fission producin g various neutron multiplici- 
ties ( J. TerreHlll957| ). we estimated a very conser- 
vative upper limit of 0.79 (10) fission bursts due 
to 238 U for Phase I (Phase II). This limit also 
assumed a 65% detection efficiency for neutrons 
above 4.5 MeV and a 40% efficiency for detecting 
a gamma burst in coincidence with the fission. 

A very small fraction of the time, the produc- 
tion of even a single neutron via the NC reaction 
can lead to an apparent burst. If a 7 ray from 
the capture on deuterium (in Phase I) or chlo- 
rine (Phase II) produces an electron above thresh- 
old via Compton scattering and then subsequently 
photodisintegrates a deuteron, it can create a sec- 
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ond neutron. Based on the SNO Monte Carlo sim- 
ulation, this photodisintegration background has 
been estimated to cause a small background of less 
than 0.01 coincidences in Phase I and 0.43 ± 0.03 
coincidences in Phase II. 

Antineutrinos can also lead to apparent bursts, 
due to the primary positron created in the interac- 
tion and the following capture of one or two neu- 
trons. Antineutrinos from radioactive nuclei in the 
earth surrounding the detector lead to a negligi- 
ble background from this process for Phase I and 
0.5 ±0.1 bursts for Phase II. A study of all com- 
mercial nuclear reactors within 500 km of SNO 
determined that the coincidence background from 
antineutrinos from these reactors was also very 
small, 0.019 ± 0.002 coincidences for Phase I and 
1.4 ±0.3 for Phase II. 

During the SNO detector construction the 
acrylic was exposed to air containing radon. In the 
decay chain of radon is 210 Po, which can decay via 
a-emission. The a can interact with the carbon 
in the acrylic, leading to 13 C(a, n) 16 reactions 
in which the 16 will produce a e + e~ pair or a 
7-ray that can photodisintegrate deuterium. The 
estimates of (a, n) coincidences from Monte Carlo 
studies of this background are low. Coincidences 
due to the diffuse supernova background (DSNB) 
and to instrumental background events are also 
estimated to be quite low, as shown in Table [5J 

2.4. Analysis Cuts 

We developed a set of analysis cuts, beyond 
the standard cuts used by othe r SNO analy- 
ses (|B. Aharmim et q/J[201ol l2007t) . to reduce the 



level of correlated backgrounds shown in Table [2] 
and discussed in the previous section. As in previ- 
ous analyses, we utilized a fiducial volume radius 
of 550 cm, as well as a variety of instrumental 
cuts based on PMT charge and timing informa- 
tion. Our high level cuts incorporated information 
such as the isotropy of the detected light and the 
event's reconstruction quality, but we did not in- 
clude any of the cuts SNO previously designed to 
remove bursts from the data set. Instead we de- 
signed new cuts that could discriminate between 
background bursts and potential supernova bursts. 

False bursts from atmospheric reactions often 
have a high-energy primary, followed by delayed 
neutron captures. We removed these bursts by 



imposing a deadtime window following any event 
whose energy exceeded roughly 80 MeV. The 
80 MeV threshold was chosen to minimize the ac- 
ceptance loss for neutrinos from a standard super- 
nova, whose energies tend to peak near 20 MeV. 
For Phase I, in which the neutron capture time is 
~ 50 ms, we used a deadtime window of 600 ms, 
and for Phase II, in which the capture time is 
~ 5 ms, we used a deadtime window of 200 ms. 
When we applied this cut to our Monte Carlo 
simulation of a standard supernova, no genuine 
supernova bursts were removed. 

In addition to having large energies, the pri- 
maries from an atmospheric neutrino interaction 
often have multiple tracks or are heavier par- 
ticles than the electrons expected either from 
charged current supernova interactions or from 
the Compton scattering of 7 rays released in 
neutron capture. We removed false bursts as- 
sociated with these events by tagging any event 
that was not electron-like, and removing all events 
within 600 ms (Phase I) and 200 ms (Phase II) 
afterward. Our definition of an electron-like 
event was based on the PMT hit pattern and 
Kolgmogorov-Smirnov tests on the angular dis- 
tribution of Cherenkov light of reconstructed 
eve nts, as has been done in other SNO analy- 
ses ( B. Aharmim et alll2010h . We estimated the 
acceptance of this cut by applying it to our Monte 
Carlo simulation of supernova bursts, and we 
found that for both phases of the experiment, only 
1.2% of genuine supernova bursts were removed. 
The same cut applied to simulated atmospheric 
neutrino events removed 57% of those bursts that 
pass the standard analysis cuts from Phase I and 
63% from Phase II. 

With the exception of accidental coincidences, 
almost all sources of false bursts are spatially 
correlated. Neutrons from atmospheric interac- 
tions or fission will capture near the primary 
event and near one another. To remove events 
on this basis, we use a cut developed by the 
Super-Kamiokande collaboration in thei r trigger- 
less burst analysis (|M. Ikeda et all\200l\ ). We de- 
fine Ar as a weighted mean of the distances be- 
tween the reconstructed event positions in a can- 
didate burst: 



Ar = 



v M-l v Af I 



mC 2 



(1) 
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Table 2 

Physics backgrounds that could create false bursts 



Backgrounds 


Phase I 


Phase II 


Phase I 


Phase II 


Phase I 


Phase II 




N ovcnt = 2 


Nevent = 2 


Novont = 2 


Novont = 2 


Novont = 3 


Novont = 3 




0.05s window 


0.05s window 


0.2s window 


Is window 


10s window 


10s window 


Atmospherics 


4.5 ± 0.9 


6.9 ± 1.4 


1.9 ± 0.4 


0.6 ± 0.1 


1.5 ± 0.3 


7.5 ± 1.5 


Muoti Spallation 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


Fission 238 U 


< 0.8 


< 10 


< 0.8 


< 10 


< 0.8 


< 3 


Photodisintcgration 


< 0.01 


0.4 ± 0.03 


< 0.01 


0.4 ± 0.03 


< 0.01 


< 0.4 


gco-v 


0.0 


0.5 ± 0.1 


0.0 


0.5 ± 0.1 


0.0 


< 0.5 


rcactor-i/ 


0.02 ±0.002 


1.4 ± 0.3 


0.02 ±0.002 


1.4 ± 0.3 


< 0.02 


< 1.4 


(a,ne"*"e~) 


0.02 ±0.10 


0.07 ±0.07 


0.02 ±0.10 


0.07 ±0.07 


< 0.02 


< 0.07 


(n,2n) 


< 0.02 


< 0.07 


< 0.02 


< 0.07 


< 0.02 


< 0.07 


DSNB 


< 0.005 


< 0.005 


< 0.005 


< 0.005 


< 0.005 


< 0.005 


Instrumcntals 


< 0.03 


< 1 


< 0.03 


< 1 


< 0.03 


< 1 



Note. — Since the search windows have different energy thresholds, we estimate the backgrounds for each search window 
for both Phase I and Phase II. These backgrounds arc further reduced by the analysis cuts described in the text. 



where \fi — fj\ is the distance between the recon- 
structed positions of events i and j within a burst, 
M is the multiplicity of the burst, and mC 2 is 
the number of non-redundant combinations. For 
SNO, we can improve on the straightforward Al- 
eut by using the fact that most spatially correlated 
events are due to neutron captures, and therefore 
there is a direct relationship between Ar of the 
events and the time that separates them, due to 
the neutron diffusion time within the heavy water. 
We therefore examined not just the spatial sepa- 
ration of burst events but also their separation in 
time, At, which is defined similarly to Eq. [TJ Fig- 
ure [1] shows a ArAt distribution for N cvcnt = 2 
from Phase II where the two-dimensional cut is 
clearly beneficial. The figure also shows the same 
distribution for Phase I, where the longer neutron 
capture time and longer diffusion distance makes 
the two-dimensional cut less effective and hence we 
use only the cut on Ar. For both phases we applied 
the two-dimensional cut to the Advent = 3 search. 
Table [3] summarizes the fraction of expected su- 
pernova signal that survives the Ar or ArAt cut 
for each of the search windows. The functional 
form of the ArAt cut was generally box-shaped, 
though for several of the search windows the shape 
of the c ut included a round ed corner described by 
At > + a 2 a 3 — a 2 Ar. Table [3] shows the pa- 
rameter values used for the ArAt cut for each of 
the search windows. In Phase I the long diffusion 
distance makes eliminating backgrounds more dif- 
ficult, meaning that the Ar cut must be harsher 



in order to adequately reduce the probability of 
seeing a false burst. 




200 400 600 S00 10001200 200 400 600 S00 10001200 

Ar (cm) Ar (cm) 




Ar (cm) Ar (cm) 



Fig. 1. — a) Simulated background bursts in the 
1 s search window for Phase II: The region inside 
the box was removed by the ArAt cut, which elim- 
inates most of the background events; b) Simu- 
lated supernova signal in the 1 s window for Phase 
II: The ArAt cut removes very little of the sim- 
ulated signal; c) Simulated background bursts in 
the 0.2 s search window for Phase I: The region 
marked with an X was removed by the Ar cut, 
which eliminates most of the background events; 
d) Simulated supernova signal in the 0.2 s search 
window for Phase I: The Ar cut must remove a 
large fraction of the supernova signal in order to 
sufficiently eliminate background bursts. 
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Table 3 

Fraction of standard SN signal that survives the ArAt cut 





Ar Cut 


At Cut 


Oil 


a 2 


a 3 


SN Signal 




(cm) 


(ms) 


(ms) 


(ms) 2 /(cm) 


(cm) 




Phase I, 0.05s window 


626 










0.389 


Phase I, 0.2s window 


668 










0.303 


Phase I, 10s window 


600 


180 


180 


100 


376 


0.867 


Phase II, 0.05s window 


385 


40 


40 


5 


200 


0.811 


Phase II, Is window 


410 


45 








0.977 


Phase II, 10s window 


400 


35 


35 


10 


300 


0.999 



Note. — Only the Ar cut was applied to the N event = 2 searches in Phase I, while we 
applied the ArAt cut to Phase II and to the N event = 3 search in Phase I. The ArAt cuts 
were generally box-shaped with one rounded corner described by the a parameters. 



3. Results 

After estimating all of the background sources 
of false bursts that could mimic a supernova 
burst signal, we placed our analysis cuts to en- 
sure that we expected no more than 0.11 events 
in each search window. Our approach to the low- 
multiplicity search was inherently blind: we did 
not examine the real data until we had set all of 
our search parameters, and we did not change any 
of those parameters after we performed the search. 

Our ArAt cut defined an "antibox" region that 
was excluded from the supernova burst search, as 
shown in Fig[TJ Prior to opening the box and ex- 
amining the events passing the analysis cuts, we 
examined the antibox in order to confirm our es- 
timates of background bursts. We found approx- 
imately the same number of bursts there that we 
expected, as shown in Table |3 Because our back- 
ground was dominated by atmospheric neutrino 
bursts and several of our other background es- 
timates were conservative upper limits, Table |3] 
compares the number of bursts we observed out- 
side the box to the number of bursts we expected 
due to atmospheric events. 

After all the cuts were developed and tested 
on simulations, we fixed our analysis and per- 
formed our burst searches on both the Phase I 
and Phase II data. We observed no bursts in any 
of our search windows, as summarized in Table [5] 



Because our 10 s window search was optimized 
for an iV cvcnt = 3 burst, we did observe some 
N eV ent = 2 bursts in that window. Using the 
same energy threshold and analysis cuts designed 
for the Phase II 10 s N CV cnt — 3 search, we ob- 
served 14 Advent — 2 bursts, which is in keep- 
ing with our expectations from accidental coin- 
cidences. Figure [2] shows the distribution of Ar 
and At for these iV even t = 2 bursts. The distri- 
bution in Ar is approximately uniform, which is 
consistent with the hypothesis that these are acci- 
dental bursts, though it is also consistent with the 
hypothesis that these are supernova bursts. The 
At distribution supports the conclusion that these 
are accidental bursts since the events are spread 
evently throughout the time window, in contrast 
with our expectations from an actual supernova. 
The event with the lowest At separation also has a 
Ar too low to have survived any of the ArAt cuts 
for the multiplicity two searches, meaning that it 
is more likely to have been a background burst 
than a genuine supernova burst. We also observed 
two iV eV ent = 2 bursts in the AT event = 3 search 
for Phase I. These two bursts are separated by 
3.2 s and 8.0 s respectively, which puts them well 
outside of the -ZV ovon t = 2 search windows. These 
bursts are also consistent with our accidental co- 
incidence expectations. 
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1000 2000 3000 4000 5000 6000 7000 8000 9000 

At (ms) 



Fig. 2. — Ar and At distributions for Advent = 2 
bursts found in the N even t = 3 search window 
for Phase II. These bursts are consistent with ex- 
pected accidental coincidences and are not consid- 
ered candidate supernova bursts. 



4. Discussion and Conclusions 

We performed a triggerless search for low- 
multiplicity bursts in data from Phase I (D 2 
only) and Phase II (D 2 loaded with NaCl) of the 
Sudbury Neutrino Observatory, finding no candi- 
date bursts. For the period of overlap our results 
are consistent with the null signal observed by 
Super-Kamiokande. 

Figures [3E] show the sensitivity of our various 
search windows to a standard supernova for both 
phases as a function of supernova distance. Those 
figures do not include the small dispersion effects 
arising from the nonzero mass of the neutrino. For 
Phase I, which was completely overlapped by the 
run time of Super-Kamiokande, our search was 
primarily looking for non-standard supernova sig- 
nals in which v e emission would be suppressed, 
allowing SNO to detect a neutral current signal 
that Super Kamiokande might not have seen. 

At a typical distance in our galaxy, 10 kpc, we 
retain a 100% detection probablity for a standard 
core-collapse supernova. In Phase I we maintain 
a 50% detection probability for a standard super- 
nova out to 60 kpc. In Phase II we retain a 100% 
detection probability out to 30 kpc and a greater 
than 50% detection probability out to 70 kpc. 
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Fig. 3. — Probability of detecting a supernova in 
Phase I data assuming a standard Burrows model. 
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Fig. 4. — Probability of detecting a supernova 
in Phase II data assuming a standard Burrows 
model. 
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Table 4 

Low-multiplicity bursts in the background-contaminated region 



Search description 


Bursts expected 


Bursts found 




in ArAt 


in ArAt 




excluded region 


excluded region 


Phase I, 0.05s window 


4.53 ±0.91 


2 


Phase I, 0.2s window 


1.90 ±0.38 


1 


Phase I, 10s window 


1.51 ±0.30 


2 


Phase II, 0.05s window 


6.94 ± 1.39 


7 


Phase II, Is window 


0.57 ±0.11 


1 


Phase II, 10s window 


7.51 ±1.50 
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Table 5 

Low-multiplicity bursts found in the SNO dataset 



Search description 


Bursts expected 


Bursts found 




in search region 


in search region 


Phase I, 0.05s window 


< 0.11 





Phase I, 0.2s window 


< 0.11 





Phase I, 10s window 


< 0.11 





Phase II, 0.05s window 


< 0.11 





Phase II, Is window 


< 0.11 





Phase II, 10s window 


< 0.05 






Note. — We observed no candidate bursts in any of the search 
windows. 
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